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Abstract: We apply the Linked Dipole Chain model to hadronic collisions using a modi- 
fied version of the LDCMC Monte Carlo program. In particular we investigate the effects of 
multiple scatterings, which in this framework are reformulated in terms of the production 
of multiple chains of dipoles. We find, among other things, that this way of describing 
multiple scatterings is less sensitive to the treatment of soft, non-perturbative, scatterings. 
Although the results presented here are only at parton level, by comparing with corre- 
sponding results from Pythia, we are confident that our model will be able to describe 
essential features of high energy hadronic collisions, including underlying events and jet- 
pedestal effects in high p_\_ events. The model should also be applicable to estimates of 
transverse energy flow in high energy heavy ion collisions, of interest in analyses of signals 
for quark-gluon plasma formation. 
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1. Introduction 

Although the main aim for the future Large Hadron Collider at CERN is to discover and 
study the Higgs boson, and to possibly discover new physics, it is important to keep in 
mind that our understanding of the strong interactions, which make up the bulk of the cross 
section in hadronic collisions, is far from complete. The hard production of large trans- 
verse momentum jets may be described to a reasonable precision by perturbative QCD, 
but for the dominating soft and semi-hard interactions we have to rely on phenomenologi- 
cal models. There may be some gold-plated signals for new physics, where these so-called 
minimum bias interactions may be disregarded, but as soon as observables include jets, a 
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more detailed understanding of the event structure is needed. Not only is the luminosity 
so high at the LHC that any signal event will be accompanied by several overlayed min- 
imum bias collisions. In addition the inclusive cross section for parton-parton scatterings 
becomes larger than the total non-diffr active cross section already at present day energies, 
which implies that there are several parton-parton scatterings in each single hadron-hadron 
collision. Such additional scatterings can be rather soft, and can therefore not be described 
by purely perturbative calculations. The production of minijets is also of interest in high 
energy heavy ion collisions, and the analysis of signals for plasma formation. They are 
essential for the initial conditions in studies of the flow in either a quark-gluon plasma or 
a hadronic phase [Q] . 

Further studies of hadronic collisions are also interesting in its own right, in order 
to understand the nature of strong interaction and the interplay between soft and hard 
physics. Many models have been presented to describe and simulate hadronic collisions. 
The Fritiof model successfully describes hadron collisions at lower energies, up to the 
ISR range. It is based on a soft momentum transfer and gluon bremsstrahlung due to the 
separation of colour charges. These charges are in the model connected by strings of the 
same nature as in e^e~-annihilation. For higher energies hard sub-collisions were added 
The formalism did, however, not give a unique description for how to combine the soft and 
hard subprocesses, and therefore the predictive power at high energies is reduced. 

The Dual Parton Model (DPM) [|, || aims at a combined description of soft and hard 
interactions in an eikonalized form, including besides a soft and a hard pomeron also triple- 
and loop-pomeron contributions. The aim is to describe not only non-diffractive events, but 
also elastic scattering and single and double diffraction, and the model is able to reproduce 
a large set of experimental data 

At high energies, in the pp collider regime, the cross section for hard or semi-hard 
parton sub-collisions becomes large, and the properties of non-diffractive events appear 
to be dominated by these subprocesses. In this energy range it is therefore essential to 
have a good description of jets from hard parton collisions, including initial and final state 
bremsstrahlung. In ISAJET and HERWIG a description of the hard sub-collisions is 
complemented by an "underlying event" due to interactions between the beam remnants. 
In Pythia I^] this underlying event is modeled by including multiple parton interactions 
in the generation of the hard parton collision, rather than adding them as an independent 
component to the hard sub-collisions. 

Even if the parton-parton collisions are of perturbative nature, with a p± above a 
couple of GeV, there are many difficulties, which cannot be solved by perturbative methods 
alone. The parton sub-collisions are typically modeled using collinear factorization. This 
gives a divergent expression for the inclusive cross section at small p±, which has to be 
treated by some kind of cutoff. Clearly collisions with large impact parameter must be 
suppressed due to screening effects, when the impact parameter is larger than the proton 
radius. The impact parameter dependence is also expected to influence the correlations 
between sub-collisions in a single event. It is natural to assume that central collisions have 
a higher activity than more peripheral collisions and therefore, on average, a larger number 
of sub-collisions Both these effects, the p± cutoff and the multiplicity correlations. 
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have to be adjusted to experimental data, which makes it more difficult to make reliable 
extrapolations to higher energies. 

In collinear factorization the hard partonic cross section is convoluted with parton 
densities, which are evolved using DGLAP [0, 0, g evolution. The parton densities 
are fitted to data on e.g. F2 in DIS and high E±_ jets in hadronic collisions. For moderate 
p_L at very large cm. energies, y^, the parton evolution contains large logarithms of s/p\ 
which need to be resummed to all orders. For these pj_ values it would therefore be more 
appropriate to use BFKL evolution ||l^, 16 ] , which gives the correct description, at least for 
asymptotic energies and fixed coupling. The BFKL evolution is related to A;j_-factorization, 
and in ref. |17| it is argued that this approach gives a dynamical suppression for small p_\_, 
which reduces the infrared sensitivity. This is mainly a consequence of a suppression of 
the partonic cross section when the virtualities of the colliding partons are larger than the 
transverse momentum exchanged in the hard collision. 

The study in ref. [0] is based on the Linked Dipole Chain (LDC) model [18, 19 1 for 
DIS. This model is based on the observation that the dominant features of the parton 
evolution is determined by a subset of emitted gluons, which are ordered in both positive 
and negative light-cone components. These primary gluons make up a linked chain of colour 
dipoles, and the sum over all such possible chains, convoluted with non-perturbative input 
parton densities, will give the proton structure function. (For exclusive properties of the 
final states also final state radiation has to be included.) The LDC model is completely 
forward-backward symmetric - i.e. it does not matter if the evolution is performed from the 
proton towards the virtual photon or vice versa. This means that LDC is also very suited 
for a description of hadronic collisions, where the sum over all possible chains, convoluted 
with input parton densities from both hadrons, will give the inclusive partonic cross section. 
The fact that transverse momenta can go up and down in a single parton chain implies 
that such a chain can correspond to more than one hard parton-parton sub-collision. Thus 
one type of correlations between the sub-collisions is automatically included in the LDC 
formalism. However, even so the per chain cross section exceeds the total non-difi^ractive 
cross section. Therefore there will in general be several dipole chains stretched between 
two colliding hadrons, and it is also in this approach necessary to make assumptions about 
the impact parameter dependence, now for the distribution in the number of chains in one 
event. 

In this paper we will investigate the effects of such multiple chains in more detail, 
using the LDC event generator LDCMC [pc| ], modified to generate hadronic collisions. The 
investigation will still be on a somewhat qualitative level since effects of colour connection 
between the chains, and of hadronization are not included. Therefore we will here compare 
our results with the Pythia generator on parton level rather than with experimental data. 

A very essential result from these MC studies is the fact that the number of chains 
per event is stable against variations of the p± cutoff in the perturbative parton evolution. 
When fitting experimental data on F2, the input soft gluon distribution is adjusted, and 
corresponds to a scale given by the cutoff. The soft input gluons can also form direct 
chains between the colliding hadrons, and these chains are also allowed to emit final state 
bremsstrahlung below the cutoff. If the soft cutoff is increased there will be fewer hard 
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chains, which is compensated by a more singular input gluon distribution. As we will 
show below, when the cutoff is increased, the decrease in the number of hard chains is 
compensated by a larger number of soft chains (which are now also allowed to radiate 
more). Thus the total number of chains in pp collisions can be fully determined by data 
on F2, without any additional adjustable free parameter. 

The outline of this paper is as follows. In section |2| we give some further comments 
on parton evolution in DIS, including a brief introduction to the LDC model. In section ^ 
we discuss the application of the LDC model for hadronic collisions. Then in section |^ we 
discuss some general aspects of multiple scatterings, in particular how they are included in 
the Pythia generator and how we treat them in the present analysis. The technical details 
of the Monte-Carlo implementation are discussed in section |5|. In section |6| we present some 
preliminary results, followed by conclusions in section |^. 



2. Small-x physics in DIS 



lepton 




Qn+l 



proton 



Figure 1: A fan diagram for a DIS 
event. The quasi-real partons from the 
initial-state radiation are denoted qi, and 
the virtual propagators fc^. The dashed 
lines denote final-state radiation. 



Effects of BFKL evolution have been studied 
extensively at HERA, where BFKL predicted a 
strong power-like increase of F2 for small x. It 
turns out, however, that BFKL suffers from huge 
next-to-leading order corrections |^l|] and it has 
been very difficult to obtain reliable predictions, 
especially when it comes to non-inclusive observ- 
ables involving final-state hadrons. More reli- 
able predictions for final-state observables can be 
obtained using CCFM evolution ||2^, which 
agrees with BFKL in the asymptotic limit, but 
also approximates DGLAP for large x. The cen- 
tral feature of the CCFM model is the angular 
ordering of the parton chains contributing to the 
proton structure function, which is a consequence 
of colour coherence or soft gluon interaction. The 
parton chains are ordered in energy (or positive 
light-cone momentum g+) and in angle (or rapidity 
y = ^ ln(g+/(7_)). (We use qi and fcj to denote the 
real emitted partons and the virtual links respec- 
tively, as indicated in figure ||.) Other kinemati- 
cally allowed emissions (symbolized by the dashed 



lines in figure ^ are treated as final-state emissions. 



2.1 CCFM Evolution 

Like BFKL evolution, the CCFM model is based on the /c_L-factorization formalism, and 
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the unintegrated distribution function is given by 



/n ,2 / -I 1 \ 

Yladz,^ (-Ar,e{zi,kl„q,) + A5 X (2.1) 

6{x - Uzi)0{qi - Qi-iZi-i)S{k'i - kljd{q - q^Zn) 

where the non-eikonal and the Sudakov form factors are given by 

/ - , 1 , ^ I \ . f - f dq'^ii dz ^ \ , , 

Ane = exp ( -aln - In — ] ; As = exp ( -a / -border I (2.2) 



z zq-^J \ J qii 1 

Here a = Sos/vr, the sphtting parameter z is defined as Zi = k^^i/k+^i-i, and the kinematic 
range allowed by the ordering constraints is specified by @ order- The quantity q^ is defined 
by Qi = 'l±i/{^ — Zi), which implies that the angular ordering condition is satisfied by the 
constraint q^ > q^^iZi^i, accounted for by a 0-function in eq. (|2.1| ). The CCFM model 
has been developed assuming purely gluonic chains, and only the singular terms in the 
splitting functions, proportional to 1/z and 1/(1 — z), are included in a fully consistent 
way. The 1/z pole is most important in the BFKL region for small x, while the 1/(1 — z) 
pole and the Sudakov form factor are essential for larger x and large k±. We also note 
that the distribution function, Q, depends on two separate scales. Besides the transverse 
momentum, k±, of the interacting gluon, it also depends on q, which determines an angle 
beyond which there is no (quasi-) real parton in the chain. 



A Monte Carlo implementation of CCFM, CASCADE |24|, reproduces a wide range of 
final-state observables at HERA, although there are still some unsolved questions. Since 
CCFM does not include the non-singular terms in the gluon splitting function, DGLAP is 
not fully reproduced in the relevant limit, and including non-singular terms degrades the 
reproduction of HERA data, especially for observables which are supposed to be sensitive 
to BFKL effects, such as forward jet cross sections |2[ 

2.2 The Linked Dipole Chain Model 



The Linked Dipole Chain (LDC) model [|l^, 19| is a reformulation and generalization of 



CCFM. It is based on the observation that the dominant features of the parton evolution is 
determined by a subset of emitted gluons, which are ordered in both positive and negative 
light-cone components, and also satisfies the relation 

q±i > mm{k±i,k±^i^i). (2.3) 

Li LDC this subset of "primary" gluons forms a chain of initial state bremsstrahlung 
(ISB), and all other emissions are treated as final state bremsstrahlung (FSB). 

This redefinition of the ISB-FSB separation implies that one single chain in the LDC 



model corresponds to a set of CCFM chains. As was shown in ref. [18|, when one considers 
the contributions from all chains in this set, with their corresponding non-eikonal form 
factors, they just add up to one. Thus, the non-eikonal form factors do not appear explicitly 
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in LDC, and the gluon distribution function is given by 



g{x,kl) ~ I ^c,dzi-^P{z,)l^s X (2.4) 



n 



Here P{z) is the fuh gluonic sphtting function, and the associated Sudakov form factor, 
defined by 

As = exp (^-a j ^zdzP{z)eorder^ (2.5) 

where again Qorder specifies the region of phase space allowed by the ordering constraints. 
The result of this reformulation is a simpler form for the unintegrated distribution functions, 
which essentially depends only on a single scale, /c^. As many of the gluons which make 
up the initial-radiation chain in the CCFM model are treated as final state radiation in the 
LDC formalism, most of the problem of angular ordering is postponed to the treatment of 
the final state radiation. The ordering of the LDC evolution in both q^ and q_, and the fact 
that the 1/z pole in the splitting function is here not associated with an non-eikonal form 
factor, imply that this contribution is completely forward-backward symmetric. Therefore 
the chain in figure || can be thought of either as an evolution from the proton towards the 
photon or, equivalently, from the photon towards the proton end. The 1/(1 — z) pole in 
P{z) and the Sudakov form factor are most essential in the DGLAP region, where chains 
with ordered q± values dominate and large z-values are important. In the LDC model the 
expression in eq. ( p.4p is symmetrized so that steps downwards in transverse momentum 
with large values of the backwards splitting parameter z-i = k-i/k-^i+i are weighted by 
the splitting function P{z-i) and a correspondingly defined Sudakov form factor. The 
absence of the non-eikonal form factor also implies that quark links and non-singular terms 
in the splitting function can be included in a natural and straightforward way. 

The set of "primari/" gluons in eq.( |2.4| ) make up a chain of linked colour dipoles, 
and the sum over all such possible chains, convoluted with non-perturbative input parton 
densities, will give the proton structure function. Besides this inclusive description of the 
events, the result in eq. ( |2.4| ) can also be interpreted as the production probability for 
an exclusive final state. We note, however, that due to the different separation between 
ISB and FSB in LDC and CCFM or other approaches, it is important to include final 
state emissions before comparing results for exclusive final states with results from other 
formalisms or with experimental data. 

The LDC model is also implemented in a MC event generator, LDCMC. When the soft 
input parton densities are adjusted to fit data on F2, the model is able to reproduce a 



wide range of HERA data |20|. The resulting integrated gluon distribution also agrees well 



with global fits by the CTEQ ^ and MRST [g^] cohaborations |2|]. However, as for the 
CCFM Monte Carlos, the data for forward jets can only be described if non-singular terms 
are omitted from the gluon splitting function. 

The fact that the LDC formalism is manifestly forward-backward symmetric implies 
that it automatically takes into account contributions from "resolved virtual photons". 
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Figure 2: (a) A chain of n initial state emissions in a hadronic collision, (b) An exchange of a 
single soft gluon without initial state emissions. 

where the chain describes evolution from the proton and photon ends towards a central 
hard parton sub-collision. This property is a consequence of the specific choice of separation 
between initial and final state emissions adopted in the LDC model, and it means that the 
formalism is also very well suited to describe hard interaction in hadron-hadron collisions, 
as will be discussed in more detail in the following section. 

3. The Linked Dipole Chain Model for Hadronic Colhsions 

As mentioned in section |2|, to get the gluon distribution function in a proton, the perturba- 
tive chains in eq. ( |2.4[ ) have to be convoluted with non-perturbative input parton densities, 
go{xo,kj_f^), for kj_ = kj_(): 

g{x,k\)= j dxodx'gix' ,k\;k\Q)gQ{xQ,k\Q)5{xQX - x) (3.1) 



As we are assuming a running coupling, which diverges for small g^, it is necessary to 
introduce a cutoff, k±_Q, for the g^j distributions in eq. Q (see e.g. H). In eq. this 
is explicitely indicated as an argument in Q. Thus (7o(2;o;^lo) assumed to describe all 
chains with transverse momenta below the cutoff k±Q. 

3.1 Hard Chains 

The fact that Q is completely forward-backward symmetric, as discussed in section ^, now 
implies that the cross section, cTc, for formation of a chain between two colliding protons 



is directly obtained if the distribution in eq. (3J) is convoluted with a parton distribution 



from the upper proton in figure |2[ Thus we find: 

= 7772- / dxQ+dx dxQ_go{xo+,k\Q)g{x,k\Q]k]_Q)go{xQ-,k]_Q)5{xQ+x xq- - k\_Q/s), 



2P 



(3.2) 
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with xo+ = A;o+ /-Pa+ and xq- = kn-/Ph-- The normahzation factor 7r^/2 originates from 
the fact that a hard scattering in the middle of the chain, with g^- ~ « j+i) is in 
eq. ( p.4D given a weight 



(3.3) 



where t ~ — ^Ij- This should be compared with the well-known cross section for gluon- 
gluon scattering, which for s » \t\ is given by 

The factor follows because Q is defined as a density in \uk\, rather than in fc^. 

3.2 Soft Chains 

It should be noted that o"c gives the cross section for a chain with at least one link with 
virtuality larger than /c^q. There will also be chains with no perturbative emissions above 
this cut. (Emissions below the cut have to be added later as final state radiation.) The cross 
section for these soft chains, i.e. the cross section for exchanging a single (non-perturbative) 
gluon with virtuality fc^Q (see figure ^), is obtained by simply convoluting the two input 
gluon densities: 



<7co(^io) = / dxo+dxQ-g{xQ+,k\Q)g{xo-,k\Q)5{xQ+XQ^ - k^^/ s). (3.5) 



^±0 

As we will see in section if the cutoff is increased, the cross section for the perturbative 
chains will decrease, but the cross section for the non-perturbative chains will increase since 
the fitted input gluon density is then more divergent for small x. This implies that the 
total number of chains in pp collisions is insensitive to the cutoff A;_lo, and therefore can be 
fully determined by data on i<2, without any additional adjustable free parameter. 

The results in eqs. ( |3.2| ) and ( |3.5| ) give together with the total non-diffractive cross 
section the average number of chains in one event. This average number is, however, not 
enough to determine the properties of exclusive final states. The sub-collisions may be 
correlated so that central collisions have a larger, and peripheral a smaller, number of 
chains. These correlations will be studied in more detail in section |^. 

3.3 Final State Bremsstrahlung 

In the LDC model, final state radiation is emitted from the dipoles in the chain in the same 
way as in the Dipole Cascade Model (DCM) p^ , |30[ for time-like parton cascades, which 
is implemented in the ARIADNE program [^]. In these time- like cascades the emissions 
are, however, limited by the k±_ of the propagators in the initial chain. In figure I, this 
corresponds to the area below the line connecting the qi points. 

The parameters in the DCM have been fitted, together with the parameters in the 
string fragmentation [^] in Pythia, to reproduce a wide range of final-state observables 
at LEP to an extraordinary precision. The main parameters in ARIADNE are Aqcd and 
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the cutoff in transverse momentum, p^cxaf which have been fitted by the DELPHI collab- 
oration to 0.22 GeV and 0.6 GeV respectively The fact that Pj_q^i typically is smaller 
than the cutoff for the ISB, k±o, obtained by fitting LDC to F2 data, means that FSB will 
be allowed also below k±Q down to the hadronization cutoff, p_Lcut- This means that if k±Q 
is increased to k'j^Q, some gluon emissions in the initial state, such as qi in figure ^, will no 
longer be allowed, but will then instead be allowed as FSB. In this way the distribution of 
the final state partons should be quite insensitive to the precise value of k±Q. 

Special care must be taken when radiating from the dipoles connecting the first ISB 
gluons with the hadron remnants. This problem is especially important for the soft chains, 
where FSB will also be allowed below k±Q from a dipole between the remnants. This will 
be discussed further in section |^. 



4. Multiple scatterings in hadronic collisions 



At high energies the inclusive cross section for parton-parton sub-collisions becomes larger 
than the total non-diffractive cross section, which means that there are on average more 
than one hard collision per event. Indications for multiple hard collisions are observed both 
in hadronic collisions |1C] and in photoproduction |34]. The properties of exclusive final 
states then depend sensitively on the distribution in the number of sub-collisions in a single 
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event. At high energies two colhding hadrons are Lorentz contracted to two flat pancakes. 
When they pass through each other the different regions are causahy disconnected, and it 
is then natural to assume that the sub-collisions are independent for a fixed value of the 
impact parameter, b. This implies an eikonalized description, where the probability to have 
n sub-collisions is given by a Poissonian distribution with an average, n(6), which depends 
on 6, or alternatively by such a distribution with the value n = excluded. Models of this 
type include Pythia, HERWIG \ and the Dual Parton Model. 

4.1 Impact-Parameter Dependence 

We note that the result depends not only upon the average number of hard sub-collisions 
for fixed impact parameter, n(b), but also on the probability for an interaction, Pmtib), for 
this value of b. These quantities are defined so that the inclusive sub-collision cross section 
is given by ahard = f d'^bh{b) and the total non-diffractive cross section by the relation 
Ctoi = / d'^b Pint{b). Expressed in these quantities we have 

n{b) = n{h)/Pint{b) (4.1) 



In ref. |1C] it is found that e.g. the "pedestal effect", i.e. the fact that the background 
is larger around a high E±_ jet than in a minimum bias event, and the distribution in the 
number of jets in a single event can be well reproduced if the 6-dependence is described 
by a double Gaussian distribution, corresponding to the presence of a kind of "hard core", 
while a single Gaussian would severely underestimate the fluctuations. (By a Gaussian 
distribution we understand here a distribution which is Gaussian in the variable b, and 
consequently exponential in the density iP'.) 

In this paper we will not present a full investigation of such correlations between sub- 
collisions, but only study some qualitative features. The fact that there often are more 
than two jets in one single chain in our approach, will imply a positive correlation between 
the sub-collisions in addition to the one originating from the impact parameter dependence. 
Therefore the fluctuations in the number of chains due to the impact parameter depen- 
dence should be weaker than the fluctuations in the number of sub-collisions in the models 
mentioned above. To illustrate the effect of these two contributions to the correlations, we 
will study two different cases: 

• Completely independent chains, where the number of chains is determined by a Pois- 
sonian distribution 

- Tl 

Pn = e-"— (4.2) 
n\ 

where n is the 6-independent average number of chains in an event. 

• An impact parameter dependence which results in an exponential distribution in n(5) 
when weighted by the probability Pmtib)- The probability P{n') for a collision with 
a centrality corresponding to n{h) = n' is thus assumed to satisfy 

P{n) = — [ dHPmt{b)6{n{b) - n') = ne""'/" (4.3) 



^ Using the add-on JiMMY program |34 
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Figure 4: The distribution in number of scatterings for non-difFractive events in Pythia. The full 



histogram is PYT4 and the dotted histogram is PYTi (see subsection 6.1 for notation). The full 
and dotted curves are geometrical and Poissonian distributions respectively with the same average 
as the corresponding histograms. 



Here n is the average n'-value. Assuming a Poissonian distribution in n for fixed 
value of b, i.e. for fixed value of n', we then find a geometric distribution with the 
same average, n, for the number of chains, n, in a random event. For the distribution 
in the number of chains we then find the geometric distribution in n, with the same 
average, n: 

P„ = / dn'Pin') . e-"'^ = ^(^)" (4.4) 
J n] 1 + nl + n 

In figure ^ we compare the simple geometrical distribution with the distribution in 
number of scatterings in PYTHIA with a double Gaussian 6-distribution, and we find that 
the geometrical distribution actually is extremely similar. Also shown is PYTHIA with 
independent scatterings and a Poissonian distribution for comparison. 

A reason to concentrate on these two alternatives in this preliminary study, is that the 
Poissonian distribution will act as a reference for the effect of further correlations, while 
the geometric distribution can mimic the fluctuations in the favoured PYTHIA version. We 
also note that both these distributions are chosen so that they depend only on a single 
parameter, the average number of chains, n. 



5. Monte Carlo implementation 

The LDC event generator LDCMC, is written to model deeply inelastic lepton-hadron scat- 
tering or, more specifically, 7*-/i scattering. To use it for generating single dipole chains 
in hadronic collisions we lower the virtuality of the 7* to below A;_lo and use a reweighting 
procedure to obtain the correct input parton densities for the hadron and to account for 
the fact that the photon only couples to quarks, while the hadron also has a gluon density. 
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Since we expect that there wih be several chains in each hadronic colhsion, we need to 
discuss two problems concerning how to combine the generated chains: we need to consider 
the colour connections between the chains and the hadron remnants, and we need to ensure 
energy-momentum conservation, so that the chains do not take more than the available 
energy in the colliding hadrons. 

5.1 Colour Connections 

We do not expect all individual chains to be connected to the hadron remnant - this would 
give a large set of overlapping strings producing too many hadrons in the full rapidity 
region. Instead we expect that ISB partons in the ends of each chain may be connected 
to corresponding partons in other chains. One could also imagine situations where two 
chains are different branches on the same cascading tree. This possibility will, however, be 
neglected in the present analysis. To be less sensitive to the colour connections we will here 
only consider the properties of the partonic final state. However, although less sensitive, 
also the final state radiation will be affected to some extent, since the emitting dipoles 
depend on these connections. 

Clearly, the question of colour connections is a very complicated one, and we will defer 
a detailed investigation to a later publication. Here we will simply assume that the FSB can 
be added to each chain independently, also allowing the dipoles between the ISB partons 
in the ends of the chain and the hadron remnants to radiate. This will certainly somewhat 
overestimate the amount of soft partons in the events and the results presented below must 
therefore be considered somewhat preliminary. 

The FSB emissions from the remnant dipoles also pose another problem. The emission 
probability from a dipole is basically flat in Inp^ and rapidity. This means that the 
emitted gluon can easily take a large fraction of the remnant momentum. The standard 
way of dealing with this in ARIADNE is the Soft Radiation Model |^5|, but here we will 
use a different method, to be as insensitive as possible to the cutoff, k±Q. Looking at the 
initial-state emission of the gluon qi in figure it will also have a basically fiat rapidity 
distribution, except there will be a (1 — xo+)^ suppression from the input gluon density 
(see eq. ( [O] ) and table |l] in section |6.lD . If we now increase the cutoff to k'j_Q , gi should 
instead be allowed as a final state emission with approximately the same distribution. We 
have therefore chosen a procedure where for each gluon emission generated according to 
the Dipole Cascade Model in a remnant dipole, a corresponding new x'q_^ is calculated from 
the energy loss of the remnant, and the emission is only accepted with a probability 



The correct procedure would, of course, be to use the actual ratio of the relevant input 
parton densities. We have, however, checked that the results are insensitive to the details 
of the suppression (e.g. changing the exponent from 4 to 3) and will therefore in this paper 
be satisfied with this simplified treatment. 
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5.2 Energy Momentum Conservation 

The treatment of emissions from the remnant dipole is also related to the question of 
energy-momentum conservation when combining several chains. If each chain takes a 
large fraction of the energy of the incoming hadron, there will not be enough energy to have 
many chains. To ensure energy-momentum conservation we have adopted the following 
procedure when combining chains into single events. 

1. First we choose n, the number of chains to be included in the event, from either the 
Poissonian distribution in eq. ( |4.2| ) or the geometrical one in eq. ([4.4|). 

2. Then we choose n chains, either hard chains with both ISB and FSB partons or soft 
ones with only FSB, according to the values of Uc and aco- 

3. The selected chains are then ordered in the transverse momentum of the hardest 
parton. This procedure is analogous to the treatment in the Pythia model. It 
implies that the hardest chain is always accepted in the following veto procedure, 
which is necessary to ensure that the spectrum of high Ej_ jets is correctly described. 

4. Each chain is now tested in order of decreasing hardness. The sum of the positive and 
negative light-cone momentum already used up by the previously accepted chains is 
calculated and from that, new effective values of x'q_^ of the current chain are obtained. 
Chains which would take more than the available energy are vetoed. Kinematically 
allowed chains are accepted with a probability 



to approximate the change in the input parton densities after some partons have 



be more correct to here use the ratio of the actual input densities used. (Again we 
have checked that our results are not very sensitive if the exponent is changed from 



5. We are now left with riacc < n accepted chains which are combined into a partonic 
event. 

The procedure presented here is at present rather cumbersome and slow and some 
results presented here have rather large statistical fluctuations. The absence of hadroniza- 
tion also makes the current version of the program unsuited for direct comparison with 
experiment. A new publicly available version of LDCMC capable of generating hadronic 
collision events will therefore be postponed until the issues of colour connections have been 
addressed in more detail. 

6. Results 

In this section we will study features of jet and minijet production, which are sensitive to 
the properties of multiple interactions. We will look at inclusive jet production, minijet 
distributions, and the pedestal effect and other properties of underlying events. 




(5.2) 



already been extracted. As for the treatment of remnant dipoles in eq. it may 



4 to 3.) 
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1.17 
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0.55 
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4 
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Table 1: The result of the fit of the parameters for the input parton densities to data from HI 
|§, ZEUS 1^, NMC |§ and E665 [|9) in the region x < 0.3, > 1.5 GeV^. The last columns 
gives the over the number of fitted data points, respectively. Parameters in bold face have not 
been fitted. 



As example we have chosen pp collisions at the Tevatron energy, 1.8 TeV in the cms. To 
be specific we define jets with a simple cone-algorithm^ using an ideahzed calorimeter with 
cells distributed uniformly in pseudo-rapidity and azimuth angle with the size 5r] x 6(p ^ 
0.2 X 0.2. The transverse energy of all generated partons are collected in the corresponding 
cells, and any cell with a summed E± above 1 GeV is tried as a jet-initiator. If the summed 
E± within a radius of 0.7 in the r]-(p plane around such an initiator is above 3 GeV, this 
E± is assigned to a jet in the E± weighted direction of the cells included. Only jets with 
|?7| < 2.5 are included in the following. 

As discussed above we will in this preliminary study not include hadronization effects, 
which are sensitive to how the chains are colour connected. We will therefore only com- 
pare with parton distributions generated by the Pythia generator, and make no direct 
comparisons with experimental data. 

6.1 Fitting procedure 

The parameters used in the MC simulations are fitted to DIS data. To get quantitative 
results from LDC we must convolute the dipole chains with input parton densities in the 
proton. These densities are not a priori known, but must be parameterized in some way 
and fitted to data. The form of these densities are taken to be 

xn{x,klQ) = Aix''^{l-x)''\ (6.1) 

where i = dy,u.u,g and s for the d- valence, u- valence, gluon and sea-quark densities respec- 
tively (where the sea fiavour densities are assumed to be = = 2/5). The parameters 
Ai,ai,bi and the perturbative cutoff, k±o, are then fitted to reproduce the measured data 
on F2. There are some relations between the parameters which are fixed by sum rules. 
Thus and A^^ are determined by flavour conservation and Ag is fixed by momentum 
conservation. The fits to F2 do not constrain the remaining parameters very strongly, so 
we have fixed the b parameters to 3 in the valence densities and to 4 in the sea and gluon 
densities. The best fit is the one called standard in reference |28] and fit-1 in table ||. The 



fitted value of k±o is 0.99 GeV, but since we would like to check the dependence on the 
cutoff, we also present additional fits in table |^ where k±Q has been fixed to higher and 
lower values. 



PYCELL in Pythia 
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In the following we will look at a couple of options for LDC and compare with different 
cases in Pythia. (For the alternatives with a single chain (LDCq and LDCq) or a single 
sub-collision (PYTq) we can only study inclusive cross sections, and the reduction from 
energy conservation due to overlayed chains or collisions is not included.) 

LDCq: Single LDC chains with the default (fit-1) tuning. 

LDCp: Uncorrelated multiple LDC chains (fit-1) according to the Poissonian distribution. 

LDCg: Multiple LDC chains (fit-1) according to the geometrical distribution. 

LDC^: As LDCo but with the A:_lo = 1-3 GeV tuning (fit-3) (similarly for LDC^). 

PYTq: Pythia without multiple scatterings selecting low-p± events without hadroniza- 
tion. Everything else default. 

PYTi: As PYTq but with uncorrelated multiple scatterings (with default values^). 

PYT4: As PYTq but with multiple scatterings according to a double Gaussian distribu- 
tion (with default values). 

PYT^: As PYT4 but with a higher soft regularization parameter (2.5 GeV instead of 
2.1 GeV). 

Since PYT4 is known to give a fair description of experimental data, we will in the 
following use it as a reasonable approximation to nature. We should keep in mind, however, 
that the results presented here are on parton level only, and that many of the observables 
studied have not been measured experimentally, so deviations from PYT4 do not imply 
that our model is wrong. We will also mostly compare LDCp with PYTi, which both are 
based on the assumption of uncorrelated scatterings; and LDCg with PYT4, since we have 
shown in figure |^ that PYT4 is well approximated by a geometrical distribution. Note, 
however, that the procedure for energy conservation in LDC severely limits the possibility 
to have many chains, and LDCg deviates substantially from a geometrical distribution, 
as can be seen in figure ^. The effects of energy conservation in PYTHIA is much smaller 
since the ISB is modeled with a DGLAP parton shower with strongly ordered virtualities 
which means that only very rarely are partons emitted close to the hadron remnants. 

6.2 Chain multiplicity 

First we want to study the cross section for the formation of a chain in a hadronic collision. 



which is given by eqs. (|3.2D and (3^). For fit-1 at a cm. energy of 1.8 TeV, we obtain 
for perturbative chains a cross section of Cc = 315 mb, which is much larger than the 
non-diffractive cross section at the Tevatron^, find ~ 39 mb. In figure |6| we see how this 
cross section per chain varies with the soft cutoff used in the fit to DIS F2 data. The result 
is an almost linear dependence on the cutoff k±Q. (Note that the input parton densities 
have been re-fitted for each value of A;_lo.) In figure |^ we also show the cross section for 



^The indices for the PYTHIA options refer to the setting of the switch, MSTP(82), selecting the multiple 
interaction model. 

*We take the value as estimated by a simple parameterization in PYTHIA based on ref. m%. 
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Figure 5: The distribution in number of cliains per event in LDC. The full histogram is LDCg 
and the dotted histogram is LDCp. The full and dotted curves are geometrical and Poissonian 
distributions respectively with the same average as the corresponding histograms. 
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Figure 6: The cross section per chain in the LDC model as a function of the cutoff, fc^Q. The 
dashed line is the cross section for chains with at least on emission above the cutoff, the dotted line 
is for soft chains without emissions above the cutoff, and the full line is the sum of the two. Note 
that the input parton densities have been re-fitted for each value of k±o. 



formation of soft chains, represented by the exchange of a single (non-perturbative) gluon 
with virtuality /c^q (cf. eq. ( |3.5D ). We see that if the cutoff is increased, the cross section 
for the perturbative chains will decrease, but the cross section for the non-perturbative 
chains will increase since the fitted input gluon density is then more divergent for small x. 
From figure ^ we see that the sum of the two cross sections is almost independent of the 
cutoff k±Q. Thus we conclude that the total number of chains in pp collisions can be fully 
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(a) (b) 




Figure 7: The differential i?_Lmax cross section. In (a) i?j_max is the largest E± jet in a chain or 
hard scattering. Full line is LDCq, dashed is LDCq and the dotted line is PYTq. In (b) S^max is 
the largest E± in an event. Full line is LDCg, long-dashed is LDCp, short-dashed is PYT4 and 
the dotted line is PYTi 

determined by data on F2, without any additional adjustable free parameter. 

Taking the sum of dc and aco (547 mb for k±Q = 0.99 GeV) we would then expect on 
the average around 14 chains per collision. However, after the veto procedure to conserve 
energy and momentum we are left with an average of 7 and 6 chains per event for LDC p 
and LDCg respectively. 

6.3 Inclusive jet production 

In figure l^a we look at the inclusive differential cross section for the highest E± jet in a single 
LDC chain and compare it to the corresponding cross section in PythiA without multiple 
interactions (PYTq). We find that the LDCMC result is close to PYTHIA. Comparing 
LDCq and LDCq we find a weak dependence on the cutoff. The reason for this may 
be that the input densities in LDC has been fitted to F2 data which are only indirectly 
sensitive to the gluon contribution. Indeed, in LDCq the gluons carry a smaller fraction of 
the proton momentum which could account for the differences in figure |^. In the future, 
when LDCMC can be compared directly with high-i?_L jet production, it should be possible 
to make a more global fit of the input parton densities, which would better constrain the 
gluon distribution. Another reason for the dependence could be the differences in FSB 
caused by our simplified treatment of the colour connections of the chains. 

In figure |^ we show the differential -E_Lmax cross with multiple interactions added with 
Poissonian (LDCp) and geometrical (LDCg) distributions for the number of chains in 
one event. For comparison we also show PYTi and PYT4. We find that the LDC curves 
more or less lie in between the PYTHIA ones, which indicates that the procedure chosen to 
combine chains is a reasonable one. We note, however, a difference for small -E_Lmax values 
where LDCg lies significantly above PYT4. 
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Figure 8: The rapidity distribution of minijets (per event with at least one minijet) In (a) the 
full line is LDCo, dashed is LDCq and the dotted line is PYTq. In (b) the full line is LDCg, 
long-dashed is LDCp, short-dashed is PYT4, dash-dotted is PYT3 and the dotted line is PYT4. 
(The different LDC curves have been obtained using the same sample of chains which is why the 
statistical fluctuations are correlated.) 



6.4 Minijet distributions 

As another check that our procedure is reasonable, we look at the rapidity distribution 
of minijets. In figure ^a we show the results for LDCq, LDCq and PYTq. We find 
that the curves agree well, and the cutoff dependence in LDC is negligible. The slight 
asymmetry in the LDC distributions is due to the fact that the recoil strategy for the final 
state emissions from the parton closest to the photon in the LDCMC for DIS is not fully 
consistently modified in this preliminary application to hadronic collisions. 

In figure ^b we show the results for the multiple interaction models. Here the LDC 
curves are above the corresponding Pythia ones. Although there are no obvious reasons 
for the curves to agree, we note that the soft FSB in LDC may be somewhat overestimated, 
especially for the soft chains. We also note that the statistics for the LDC curves is poor 
due to the present inefficient simulation procedure. In addition, the statistical fluctuations 
in LDCp and LDCg are correlated since the two curves are obtained using the same set 
of generated single LDC chains. 

To investigate the fluctuations we show in figure ^a the distribution in number of jets 
in the central 5 units of rapidity for the different multiple interaction models. It is clear 
that LDCg has much larger fluctuations than has LDCp, as was expected. LDCp has 
also larger fluctuations than PYTi although both are based on uncorrelated scatterings. 
This is because in LDC there is, as mentioned above, a possibility to have several hard 
scatterings in each chain, which increases the correlations. LDCg does not have as large 
fluctuations as PYT4, which also was expected since, as shown in figure |5|, the distribution 
in number of chains in LDCq is much narrower than a geometrical distribution, due to 
energy conservation. 
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Figure 9: The distribution in number of minijets in the range \r]\ < 2.5 (per event with at least 
one minijet). (a) shows all minijets, while (b) is all minijets in the minimum azimuth region (as 
defined in the text) in each event with iJj^max > 10 GeV. In both cases the full line is LDCg, 
long-dashed is LDCp, short-dashed is PYT4, and the dotted line is PYTi. 

6.5 The Pedestal Effect and Underlying events 

The contribution from hard sub-collisions is in many analyses added on top of an indepen- 
dent "underlying event". In reference ||4l| a number of observables were presented which 
were shown to be especially sensitive to this background activity. The regions of azimuth 
which are transverse to the highest E± jet in an event given by 60° < (j) < 120°, should be 
less affected by the primary trigger scatterings, which mainly populate the collinear and 
opposite azimuth region from the highest E± jet with hadrons. There are, however, also 
effects from bremsstrahlung from the primary interaction in the transverse regions, but the 
authors argue that these typically only affect one of these two regions. Thus, by looking 
at the activity in the transverse region which has the smallest activity, one is especially 
sensitive to effects of the "underlying event" and secondary scatterings. We will call this 
region the minimum azimuth region, and in |9|b we show the distribution in number of jets 
in this region in events with a trigger jet of -E^max > 10 GeV. Clearly, the general trend 
found for the overall distribution in number of jets is present also here. 

One of the most important features of the underlying event is the so-called jet pedestal 
effect, i.e. the fact that the underlying event activity is not independent of a large E± jet. 
In figure ^ we show the average number of minijets as a function of the E± of the hardest 
jet, both for the overall minijet activity and for the minimum azimuth region. As expected 
all models show an increasing average number with increasing -E^max since, even in the case 
of uncorrelated scatterings, there is a higher probability to find a large E± jet in an event 
with many scatterings than in an event with few. Also, the models with most fluctuations 
gives the strongest rise, while models with fewer fluctuations saturate earlier. We note, 
however, that for this variable there is a particularly large difference between the LDCg 
and PYT4 models. 
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Figure 10: (a) The average number number ( 
largest jet in the event. The full line is LDCg, 
the dotted line is PYTi. (b) Is the same but 
region. 
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Figure 11: (a) The distribution in number of minijets in the minimum azimuth region in events 
with -B_Lmax > 10 GeV. (b) The average number of minijets per event in the minimum azimuth 
region as a function of the E±ni&x- In both plots the full line is LDCg, long-dashed is LDC^, 
short-dashed is PYT4 and the dotted hue is PYT^. 



6.6 Cutoff Dependence 

Finally we would like to check the cutoff dependence also for our multiple scattering model. 
In figure |ll| we show the distribution in number of jets and the average number of jets 
as a function of £^_Lmax for LDCg and LDCq and we find that the differences indeed 
are very small. As comparison we also show the results for PYT4 and PYT^. For the 
double-Gaussian model in Pythia, there is no direct cutoff, but rather a soft regularization 
parameter, p±Q, which forces the partonic cross sections and the running of to be finite 
at small scales by modifying them according to 

do-part (pi) C^O-partbi) f P\ ^ ^ 



as{pi)^as{pi+pio) (6.2) 
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We see in figure |Tl| that the result is very sensitive to the value of p±o, which has to be 
adjusted to fit experimental data. We also note that in Pythia the value of p±o is expected 
to increase with energy in a way which is not determined by the model, which means 
that it is difficult to make reliable predictions for future experiments at higher energies. 
Naturally we should not compare the variation in cutoff for LDC (0.99 — > 1.3 GeV) with 



the variation in the regularization parameter in PYTHIA (2.1 — > 2.5 GeV). From figure 11 



we can, however, see that the dependence in LDC is very weak, which means that in this 
model the result can be directly estimated from data from DIS. 



7. Conclusions 



In this article we present a new approach to describe high energy hadronic collisions, 
which is based on the Linked Dipole Chain model for DIS. The main advantage over earlier 
treatments is a much reduced sensitivity to an infrared cutoff, which in models based on 
collinear factorization has to be fitted to experimental data. In our model the parton- 
parton sub-collision cross section can be fully determined by data from DIS, which implies 
a more intimate relation between hadron-hadron collisions and DIS. We hope that this 
relation also implies increased possibilities to use data on hh collisions to constrain fits to 
parton distribution functions. Thus e.g. the connections between different parton chains in 
hh collisions may be related to gluon saturation effects in DIS. The insensitivity to a soft 
cutoff also makes it easier to make reliable extrapolations to higher energies. 

At high energies the properties of hadron-hadron collisions appear to be dominated 
by hard and semi-hard parton collisions. Although the parton-parton collisions are of per- 
turbative nature, various non-perturbative effects make a description of hadronic collisions 
much more complicated than e.g. DIS: 

(i) In a traditional collinear factorization scheme the inclusive sub-collision cross section 
is divergent for small p±, which has to be treated by a soft cutoff p±o. 

(ii) The inclusive sub-collision cross section is larger than the total non-diffractive cross 
section, which means that there may be many hard sub-collisions in a single event. 
This implies that the correlations between the sub-collisions are important. 

(iii) The scattered partons can be colour connected with each other and with the proton 
remnants in different ways by strings or cluster chains. 

In models based on collinear factorization the soft cutoff, p±o, has to be fitted to 
experimental data. The correlation between partonic sub-collisions is, in e.g. the DPM 
model, Pythia, and HERWIG, described in an eikonalized form by an impact parameter 
dependence. For a single hard sub-collision the colour structure can be read off from the 
perturbative matrix elements |^, but for events with many sub-collisions the result 
necessarily depends on non-perturbative effects. 

The LDC model for DIS is a reformulation and generalization of the CCFM model. 
It is based on fc^-factorization and the inclusion of non-fex-ordered parton chains, which 



- 21 - 



are expected to be important for limited and high energies. The formulation is fully 
symmetric with respect to the photon and the proton ends of the parton chain, which 
makes the model also suitable for a description of hadronic collisions. The input soft gluon 
distribution, go{x,QQ), is fitted to data from DIS, and when we then calculate the number 
of parton chains in pp collisions, we find that the result is actually fixed by the fit to F2. 
If the soft cutoff used in the fit is increased the number of hard chains is reduced, but 
this is fully compensated by an increase of soft chains described by the input distribution 
function, qq. (These soft chains are now also allowed to emit final state radiation up to the 
increased value for the A;_L-cutoff.) 

Concerning point (ii) in the list above, we note that some correlations between hard 
sub-collisions arc included automatically, as a single non-A;j_-ordcrod parton chain can con- 
tain two or more local maxima which all correspond to separate hard sub-collisions. Besides 
this correlation due to sub-collisions in a single chain, we expect that the number of chains 
is larger in central and smaller in peripheral collisions. As in the models mentioned above, 
we want to describe this by an impact parameter dependence, and in this paper we have 
studied two different distributions, one purely stochastic and one similar to the double 
Gaussian distribution used in the Pythia model. 

The third problem, concerning colour connections, is very important for the final dis- 
tribution of hadrons, but less important for a determination of the partonic state before 
its hadronization. In this preliminary study of the model we focus on features of jet and 
minijet production, such as the inclusive jet production, minijet distributions and the jet 
pedestal effect, which are less sensitive to hadronization effects. We will postpone a more 
detailed analysis of the colour connection structures to future work. In this study we there- 
fore do not include any hadronization effects, and as a consequence we are not comparing 
our results with experimental data, but are satisfied by a comparison with parton level 
results from the PYTHIA generator. 

For the numerical estimates we are using a modified version of the LDCMC event genera- 
tor program, and as an example we study collisions at s/s = 1.8 TeV. Prom the numerical 
results we realize that it is very essential to consider the effects of energy-momentum con- 
servation. This reduces the possibility to have many chains in a single event, and it also 
constrains the final state radiation in the regions close to the proton remnants. 

On a qualitative level the results from our two versions agree well with the ones ob- 
tained from Pythia. We find that our version with uncorrected chains gives larger fluc- 
tuations than the corresponding model with uncorrected scatterings in Pythia, which is 
expected since the dipole chains already include some correlations between multiple scatter- 
ings in the same chain. Including the correlations from our geometrical chain distribution 
increases the fluctuations, but not fully to the level of the model based on a double Gaussian 
impact parameter dependence in PYTHIA. This is because the suppression of large chain- 
multiplicities due to the energy conservation procedure, is larger than the corresponding 
suppression in PYTHIA. 

As discussed above, there are still unsolved problems in the LDC model for hadronic 
collisions. This concerns in particular the modeling of colour connections between the 
chains and the hadron remnants in a single event, which needs to be analysed in much 
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more detail. In turn, this will affect the way final state bremsstrahlung and hadronization 
are handled, and also the impact of energy conservation. In particular, the simplified 
model presented here probably overestimates the soft final state emissions, since the final 
state radiation is allowed from each chain independently. With proper treatment of colour 
connections we expect a reduction of such soft emissions. 

Furthermore, the implementation of the model in the LDCMC event generator needs to 
be improved substantially, as the current procedure, where 7*-p collision events need to be 
reweighted and combined by hand into p-p events, is clearly too cumbersome. 

Despite the shortcomings of the simplified model presented here, we are confident that 
the problems can be solved and that the final model will be able to describe essential 
features of high energy hadronic collisions, including underlying events and jet-pedestal 
effects in high E± events. 
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